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Abstract 



fH ' We discuss relativistic dynamics in a random electromagnetic field 

which can be considered as a high temperature limit of the quantum elec- 
tromagnetic field in a heat bath (cavity) moving with a uniform velocity w. 
We derive a diffusion approximation for particle's dynamics generalizing 
the diffusion of Schay and Dudley. It is shown that Jiittner distribution 
is the equilibrium state of the diffusion. 

>" 

Q ■ 1 Introduction 

o: 

CNJ . The dynamics of a test (tracer) particle in a medium of some other particles 

.,— b ' is often approximated by diffusion [1] . It is an old result [2] [3] that the Vlasov 

f^ . equation describing a tracer particle moving in a random electromagnetic field 

Cn ' produced by a chaotic motion of other particles can be approximated by a 

diffusion equation. The diffusion approximation finds applications in heavy 

ion coUisions [4] [5], plasma physics and in astrophysics [2] [3] [6] [7] . A detailed 

discussion of the diffusion approximation can be found in [8][9] [10][11] [12]. 

The mathematical theory of a relativistic diffusion is not so well developed 

C^ . as the non-relativistic one (see the reviews in [13] [14]). A relativistic diffusion 

in the configuration space does not exist [15] [16]. If the diffusing particle has a 
fixed mass then an analog of the Kramers diffusion (on the phase space) without 
friction is uniquely determined. The generator of the diffusion is defined by the 
Laplace-Beltrami operator on the mass shell as has been shown by Schay [17] 
and Dudley [18]. In our earlier papers [19] [20] we have studied a friction term in 
the framework of the relativistic Brownian motion of Schay,Dudley, Franchi and 
Le Jan [21]. The friction term has been defined as a drift uniquely determined by 
the equilibrium state and by the requirement of the reversibility of the process 
(in physical terms by the detailed balance). Such a drift term has been obtained 



by other means in [22] [23] . Without a friction the relatvistic Brownian motion 
is unphysical,e.g., the particle energy grows exponentiahy fast [20]. 

FoUowing a development of the theory of the non-relativistic Brownian mo- 
tion we may expect that relativistic deterministic motion, which is chaotic in 
certain sense, can also be approximated by a relativistic diffusion. In our ear- 
lier paper [24] we have discussed the dynamics of a relativistic particle in a 
Poincare invariant random electromagnetic field. We have shown that the rel- 
ativistic diffusion of Schay and Dudley can be considered as an approximation 
of the proper time dynamics in this field. However, such a Poincare invariant 
electromagnetic field is rather unphysical. In this paper we discuss particle's dy- 
namics in a thermal (random) electromagnetic field with a distinguished frame 
of reference. We can think of the Universe filled with the CMB radiation as an 
example of a Poincare covariant system with a heat bath [25]. When describ- 
ing a random system it is more appropriate to apply a statistical (Liouville) 
description. The evolution in the proper time is determined by the expectation 
value of an exponential of the Liouville operator. The expectation value is de- 
fined by the probability distribution of the electromagnetic field. The two-point 
function is expressed by its Fourier transform G. We do not need to specify 
the thermal distribution G (it can be the Planck distribution). We consider an 
observer moving in an incrtial frame with the velocity w (this is equivalent to 
moving the heat bath and keeping the observer at rest). For such an observer 
the thermal distribution is a function of the momentum k of the electromagnetic 
modes and w [25] [26] . We restrict ourselves to the expansion of the exponen- 
tial of the Liouville operator till the second order in the electromagnetic field 
strength (the square of the Liouville operator). Assuming that the correlation 
functions of the electromagnetic modes (Fourier components) are decaying fast 
for the time larger than the correlation time Tc we can apply the Markov ap- 
proximation treating the particle's subsequent scatterings upon electromagnetic 
modes as independent. According to Kubo's theory [8]- [9] the second order dif- 
ferential operator arising in the small time expansion is the generator of the 
diffusion defining the evolution of the system also at times large in comparison 
to the correlation time Tc- We show that the form of this diffusion is determined 
uniquely by the assumptions of the Lorcntz invariance of G. This assumption 
allows to express the expectation value of the square of the Liouville operator 
unequivocally up to two constants. It follows that the generator of the diffusion 
depends only on the momenta and on the frame w. The diffusion equation in 
the proper time is a generalization of the diffusion of Schay [17] and Dudley 
[18] (corresponding to zero temperature). The particle mass is preserved by the 
dynamics. It enters in an additive way in the model. We can sum over the mass 
when computing expectation values of observables, e.g., the energy-momentum 
or entropy. We think that keeping the mass as a variable (we can fix it at the 
end if necessary) can be a useful tool in some models of diffusion. In quantum 
physics [27] the Wigner function is an analog of the classical phase space dis- 
tribution . The square of four-momentum is off the mass shell in the Wigner 



function description [27] (a matter with a continuous mass spectrum has been 
suggested recently in [28]). In an appUcation of the diffusion equation to cos- 
mology [29] [30] the mass of a galaxy can be considered as a parameter varying 
continuously. The diffusion equation, which we have obtained in a model of an 
interaction with the electromagnetic thermal radiation, could have a broader 
application to any diffusive dynamics in an environment characterized by a dis- 
tinguished reference frame (the heat bath). The diffusion in the proper time 
originated from the dynamics in a random field. In a random system the mean- 
ing of the proper time is loosing its physical meaning. We should express the 
evolution of states and observables (e.g., the energy-momentum) in the coordi- 
nate (laboratory) time. We obtain a transport equation in the laboratory time 
from the requirement that probability density should not depend on the proper 
time. There is a certain similarity between our transport equation and the one 
of refs.[31][32][33]. These authors describe a diffusion of a particle in a fluid. 
The fluid supplies a distinguished reference frame. Our form of the equation 
and the constants appearing in it result from a concrete physical model which 
has not been considered in refs.[31]-[33]. 

The plan of the paper is the following. In sees. 2-3 we repeat the formu- 
lation [24] of the kinetic equation in an electromagnetic field in order to fix 
the notation. In sec. 4 the two-point function of a random electromagnetic field 
depending additionally on the four-velocity w^ is discussed. The exponential 
of the Liouville operator is expanded till the second order in the proper time. 
Its expectation value is calculated. In sec. 5 we take the diffusion generator as 
calculated in previous sections as a basis of the relativistic diffusion theory. We 
show that Jiittner distribution is the equilibrium state of the diffusion. 

2 Kinetic equation of particles in an electromag- 
netic field 

The dynamics of a particle in an electromagnetic field is described by the equa- 
tions [34] 

V^^=P', (1) 



Vp^^^^^.p", (2) 

where /i = 0, 1, 2, 3, x denotes the coordinate and p the four-momentum . J" is 
an antisymmetric field strength of the electromagnetic field. We could express 
the mass as y/p^ without changing the meaning of eqs.(l)-(2) because it follows 
from eq.(2) that p^ — r/fj^^p'^p'^ = const as a consequence of the antisymmetry 
of Tfj,i, and 

r/^>M^P. = (3) 



with r]fj,i, = (1,-1,— 1,-1) . From cqs.(l)-(2) it follows that t is proportional 
to the proper time as 

where x'^ — ct { we set c = 1 from now on for convenience). We can eliminate 
T from eqs.(l)-(2) in favor of x" if we wish. 

In the kinetic theory of classical particles we consider a distribution /C,- of 
trajectories {x{t, y , q) , p{t, y , q)) in the phase space starting from {y,q) 

ICr{x,p; y, q) = S{x - x{t, y, q))5{p - p{t, y, q)). (4) 

We can consider more general distributions of trajectories $,- by spreading the 
initial points in the phase space with a certain probability distribution $ 

<^r{x,p) = / JCr{x,p;y,q)^{y,q)dydq. (5) 

ICt as well as $r satisfy the differential equation 



Kt with the initial condition 5{x — y)S{p — q) and $T-(x,p) with the initial 
condition ^(x,p) . We defined here a first order operator Q~^ which will be 
discussed later in more detail. In eq.(6) we assumed that po is an independent 
variable. Working with the four-momentum off the mass shell allows to preserve 
explicit Lorentz invariance. A formulation with an independent po may have 
applications to quantum phase space methods based on the Wigner function 
[27] [35] [36] . li p is on the mass shell then we skip po derivative in eq.(6). 

The density of trajectories ft in the laboratory time x" can be expressed by 

$r [27] 



^{x,p)^ j dT^,{x,p). (7) 

It satisfies the transport equation 

g+n ^ 0. (8) 

In the Liouville approach to the classical statistical mechanics we consider func- 
tions W on the phase space as observables. We may define the expectation value 
4)t-(VF) of W in the state $r as 



^^{W) = ($r,VK) ^ j dxdp^r{x,p)W{x,p)= {^,Wr) 

= J dydq^{y, q)W{x{T, y, q),p{T, y, q)). 



(9) 



The semigroup (TrVF)(y,g) = Wriy^q) — W{x{T,y,q),p{T,y,q)) is unitary in 
the scalar product (9). We have from the definition of Wr in eq.(9) 

£{T^W){y,q) = iTrgW)iy,q) = {gTrW){y,q) = ^{q^^ + T^'^q,^) 



^gwr, 

(10) 
where the generator Q of the semigroup TV = exp(TtJ) is the adjoint of Q^ 
of eq.(6) in the Hilbert space L^{dydq) ( by dy or dq we denote an integral 
over W^; vectors from R^ will be denoted by a boldface letter). In the deriva- 
tion of eq.(lO) we apply the semigroup law Tr+s = T^Tr = Tt-T^ in the form 
W{x{s,Xr,Pr),p{s,Xr,PT)) = W {x{t, Xg , Ps) , p{t, Xs, Ps)) (where {xr.Pr) IS the 
shorthand notation for the trajectory of eq.(4)) resulting from the composition 
of trajectories. Then, the equality T^-G — GTr in eq.(lO) follows by differenti- 
ation. The requirement of the r-independence of the probability distribution 
(^^T — 0) gives the kinetic equation in the laboratory time t . This is the 
same equation as the one which can be derived by an elimination of r in favor 
of t in the evolution equations (l)-(2). 

3 Random dynamics 

In this section we repeat some definitions used in our earlier paper [24] in order 
to fix the notation. We write the formula (10) for an evolution of a function of 
trajectories in a general form as 

^W^{X + Y)W, (11) 

dr 



where 



and 



X = ^p^dl (12) 



Y = 4=-F^Va^- (13) 



p 

From now on derivatives over space-time coordinates will have an index x and 
the derivatives without an index mean the derivatives over momenta. Let 



1__„.,, 1 . .„ 1 



Y{t) = CXp(-TX)y CXp(tX) - -^=T'"'{X - —=pT,p)p,{d^ + ^=Tdl). 



(14) 
Then, the solution of eq.(ll) can be expressed in the form (an analog of the 
interaction picture in quantum mechanics) 

Wr = cxp(tX)VK/, (15) 



where 

a^VK/ = Y{t)W^. (16) 

We assume that the electromagnetic field is random. We consider an expectation 
value of Wr (we could equivalently discuss $r with the same result). In general, 
we have the cumulant expansion for the expectation value (the expectation value 
over the electromagnetic field T is denoted by (.), the initial condition W is 
assumed to be independent of T) 

{W^) = exp ( /„* ds{Yis)) + i /J ds /; ds'{iAis)Ais') + A{s')Ais))) + ...)w, 

(17) 
here 

A{s) = Y{s) {Y{s)) 

If [Y{s),Y{s')] = and F is a linear function of Gaussian variables then eq.(17) 
is exact (with no higher order terms). The expansion of the dynamics (11) till 
the second order term reads 

{W^) =W + j;; ds{Yis))W + i(/; ds{Yis))fW + i((/; dsA{s)y)W + ... 

(18) 
The approach of Kubo [8] -[9] approximates the random Liouville operator in 
the exponential (17) by the t^ term in the expansion (18) (we discuss this ap- 
proximation in the next section). 

4 Diffusion in random electromagnetic fields 

In ref. [24] we have calculated the generator of the diffusion in a Poincare invari- 
ant random electromagnetic field. We have assumed that the correlation func- 
tions of the electromagnetic field depend only on the space-time coordinates. 
In this paper we discuss the correlation functions which may depend also on a 
time-like vector ui^ (without loss of generality we may assume w^w^ = 1; as will 
be seen from the computations the vector w^ plays also a role of an ultraviolet 
cutoff). We have in mind a model of a particle moving in the electromagnetic 
field of the black body radiation at finite temperature /3^^ (the heat bath). The 
Lorentz invariance of the particle system is lost if we keep the heat bath at rest. 
We can express the Lorentz invariance of the whole system (a particle interact- 
ing with the electromagnetic field and the environment) considering a heat bath 
moving with a velocity w^.In quantum field theory the correlation functions of 
the quantum electromagnetic field are usually calculated in the rest frame of 
the heat bath [37] (for explicitly Lorentz covariant calculations see [38]). Cor- 
relation functions of quantum fields are not symmetric (quantum fields do not 
commute). Hence, they do not describe a random field. However, in the high 
temperature limit the non-symmetric part can be neglected [37]. In this limit 
the correlations can be represented as correlation functions of random fields. 
We split 



where i^^j^ has zero mean value. We assume that an average (•) over T is defined 
which preserves the Poincare symmetry. This means that the two-point function 
defined by 

{F^,,{x)F^p{x')) = G^,^,,p{x - x') (20) 

is a tensor. Gim-.ap is symmetric under the exchange of indices (^[iv\ x) and 
(ap;x') and antisymmetric under the exchange fi ^ v and a ^^ p. We impose 
the Bianchi identities ( as J^ is expressed by a potential) 

d^Tp,^, + dlT^p + dlT,„ = 0. (21) 

In terms of the two-point function 

a^e"'^^'^G^.;.p = 0. (22) 

For Fourier transforms eq.(20) reads 

{Fp,{k)F„p{k')) = G^,,„p{k)5{k - k'), (23) 

where G^^-ap{k) is a tensor which must be constructed from the vectors fc^,w^ 
and the fundamental four-dimensional tensors ?7^p and e^iypo-. Hence, in general 
we could have with some scalar functions a,, (higher order tensors in k and w 
would not satisfy the antisymmetry requirements mentioned below eq.(20)) 

Gpu;(jp\'^) ^ O'lKjIp.al^u'^p ^ TJupKi^Kf, + 'Tj^pK^Ki^ '>]uiT'^p,'^p) ' Oj0^p.vap 

+a2{rip,aWukp - rip^pWuka + ri^pWfik„ - rj^aWpkp) (24) 

+a^{'np,tyWv'Wp - 'qp.pWvWa + r]upWpWa - ri„aWpWp). 

However, the Bianchi identities (22) lead to 02 = 03 = 0. The reason is that in 
eq.(22) (in the x- representation) we need in addition to d^ one more derivative 
with an index i^ or /i in order to give zero on the rhs. There are terms multiplied 
by 02 without such a derivative, hence 02 — 0. It is clear that the as terms cannot 
satisfy (22) because they are not symmetric in {p,v)- Hence, 03 ~ 0. The 
requirement of positivity of the probability measure defining the expectation 
value (20) leads to oo = ( as discussed in [24]; in the case of a quantum free 
electromagnetic field at finite temperature [37] we obtain also the representation 
(24) where only ai ^ 0). Summarizing, the two-point function of a random 
electromagnetic field satisfying the Bianchi identities must have the form 



Gf^^;apix, x') = / dkGf_i^.ap{k) exp(ifc(x - x')) 

^ J dkG[k) exp(zfi;(x x ))(T]paki^kp Tjp^pk^ka + f]^pkpka T],ycrkfj,kp), 



(25) 



Hence, 

{F^Ax)F,pix'))fj = -D^,,„pG{x - x'), (26) 
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where 

Df,^;ap = -Vt^ad^d^p + ri^pdf,d% - v^pd^d^ + v^ad^d^. (27) 

The two-point function is positive definite if and only if G{k) in eq.(25) satisfies 
the condition 

G{k) > (28) 

and G(fc)=0 if fc^ < (as we have shown in [24]). 
It follows from eq.(25) that 

iDpi,-apG){x) == rip„Typ{x) - r]fj^pT^^{x) + r]^pT^p{x) - r]^crTpp{x), (29) 

where 

Ti,„ix) = -dp^G{x) = I dkGkf,K exp{ikx). (30) 

We shall need solely 

Tpv{0) = / dkGkpki, = eWf^w^ - TTeiilpv - w^w^) (31) 

The decomposition of the rhs of eq.(31) follows from the Lorentz invariance. If 
G has the meaning of the distribution of electromagnetic modes (photons) then 
e is the energy density and tTc the pressure. The constants in eq.(31) can be 
expressed by G 

'' dkG{k){kwf, 



TT, = i I dkG{k){{kwf ~k^). 

If G{k) is vanishing for A;^ < (as is the case for a probability distribution 
of electromagnetic modes) then from eq.(31) we can derive the inequality 

e > Stt, > 0. (32) 

We could express the expectation value {A{s)A{s')), needed for a calculation 
of the diffusion, by dpd^G{—j={s — s')p) of eq.(30)calculatcd at the point (s — 

s')p ( it comes from the correlation functions of F^^'^{x j=ps,p) of eq.(14) 

calculated at s and s' with the covariance (26)). This function has more involved 
tensor decomposition than T^'' of eq.(31) because there are still some tensors 
constructed from p^ and Wp on the rhs of eq.(30). However, these terms will be 
of higher order than r^ . In the limit of small r we calculate 

A^ = lim.^o ^{{f^ A{s)dsf) - lim^^o ^ Jo ds J^ ds' {A{s)A{s') + A{s')A{s)} 

= (-r/^"T'^''(0) + r]^PT'"'{0) - r^^'PT'^'^iO) + rj''^TPP{0))^p,dp^Ppd,. 

VP VP 

(33) 



We dropped the sd^ term from eq.(14) because it does not contribute to the 
limit of smah r. If we apply the formula (31) then we obtain 

+ (e + Tr,)p-^p''d, + 2wpw''d,) ^ ^ 

where 

with 
and 



'{e+^T,)p-^iipwf-{wp)iwPp''+w''p^')+p'^wPw'') = PP^C^pPP" 

(35) 



p'P^"" = 77^'^p^ - pPp" (36) 



Cc,p = 27re77o-p - (e + n^){{wp)'^p '^rj^p + w^Wp). (37) 

It should be pointed out that without the velocity w the limit in eq.(33) would 
not exist. In fact, Tp^{{)) is infinite if G depends only on k"^ . A damping factor 
in G is necessary if the integral (31) is to exist (such a damping factor is present 
in the Planck distribution (47)). 

The approach of Kubo [8] -[9] approximates the random Liouville operator 
on the rhs of eq.(18) by the r^ term. Kubo shows [8]-[9] that the kinematic r^ 
behaviour in random dynamics at times short in comparison to the correlation 
time Tc (of the electromagnetic modes) goes into the diffusive r behaviour at 
times large in comparison to the correlation time. Kubo's argument can be 
considered as a Markov approximation ( no memory of successive steps) to the 
random dynamics. 

Another interpretation of the diffusion approximation (which can be made 
rigorous [11]) interprets the small t expansion as a weak field expansion. The 
weak field can be rescaled to a large time. The rescaling can be interpreted as a 
transformation of a microscopic behaviour of a particle system into an evolution 
at a macroscopic time scale. In more detail , eq.(2) is rewritten as 

ar ^pi 

Then, rescaling the time and an expansion in 5 lead to the diffusion with the 
diffusion tensor 

a^" = r ds{FP''{x)F''P{x - s^p))p,pp. (38) 

The integral (38) is proportional to the correlation time t^- As discussed in [8]- [9] 
even if the correlations of the random fields are not decaying exponentially then 



nevertheless the behaviour of the particle dynamics in random fields described 
by the matrix a for times small in comparison to the correlation time determines 
the behaviour for times large in comparison with the correlation time. The large 
time behaviour is diffusive with the diffusion matrix a related to a in eq.(34) 

by 

a^'^^Tca^'". (39) 

The Kubo's formulas (the analogs of eqs.(38)-(39)) are widely applied in the 
transport theory of random systems. 

As will be discussed in sec. 5 in general the diffusion generated by Aw (34) has 
no equilibrium. We need a friction term in order to stabilize the behaviour at 
large time. We would like to point out that if {d§J^^„) ^ then in the expansion 
(18) there is a drift term which can lead to a friction . Let us consider in eq.(18) 
the term which is of the first order in the field. Expanding J^(x — s—^p) in s 

we obtain 

K = /;ds(r(s)) = r{J^^,{x))^p^d^^ - r^{d^J^^,)p-^p-p''d^ + o{t^). 

(40) 
Poincare invariance implies that 

(-^M-> = (41) 

because no antisymmetric tensor can be built from e'^'^"^ , rj^'^ and w^ . As- 
suming that the second term on the rhs of eq. (40) is non-zero we have from the 
Poincare invariance 

{d'^cr^tJ.u) = -^{riat^W^ - f^auW^) (42) 

with a certain constant r. Hence, 

K = rwpp-^pi.d'' - rwf.d'' = -rP'^^Wf^d,,. (43) 

with P^" = 77^'' — p^^p'^p'^ . We can get an interpretation of cq.(42) if we 
calculate an electric current of the electromagnetic field 

J, = d^T^,. (44) 

From cq.(42) 

(^^)=-|u;^ (45) 

The current (45) results from a charge — ^ moving with the velocity w^ . We can 
obtain such a current in a finite temperature quantum field theory of interacting 
electromagnetic and complex scalar fields (j) [37]. The density matrix is 

Pfj=eM-l3P''w,), (46) 



10 



where /3~^ is the temperature, P^ is the four-momentum of the quantum fields 
and w^^ describes the moving frame. Then, calculating the expectation value 
(of the electric current) we obtain in the lowest order of the expansion in the 
coupling constant 



{J^.)(i = -iTr(^exp{-pP-w,)i^d-^(f>-^d-^^)) 



= -JdkkQ^kf,{cxp{l3kw) - l)-i = - 



(47) 



w,, 



2 '-^M 



with a certain constant r. Here, {.)p denotes an expectation value (defined by 
the trace) with respect to the state (46) . 

5 Relativistic diffusion at finite temperature 

We define the diffusion generator by the lowest order ( r^) term in the expan- 
sion (18) in proper time (including the friction term coming from a constant 
current as discussed at the end of sec. 4). We incorporate Kubo's argument that 
this operator multiplied by the correlation time generates a diffusion properly 
approximating the random dynamics. In this way we are led to the following 
diffusion equation for particles with a continuous mass spectrum ( by a certain 
abuse of notation we denote by W the expression which has the meaning of (W) 
in random dynamics of sec. 3) 

j^drW = Q^W = ^p^a;^VK + ^^a^^a^VK - rP^'''w^^^W, (48) 

where r is determined by the expectation value (42). 

The diffusion in the proper time originates from the relativistic dynamics 
in the proper time. Then, the proper time has a well-defined meaning. This 
is the time measured in the frame moving with the particle. However, if the 
particle motion is random then this frame is random as well. Such a random 
frame is loosing a physical meaning. For this reason we should express the 
kinetic transport in terms of the coordinate (deterministic, laboratory) time. 
The probability distribution evolves according to the adjoint equation (6). The 
independence of the proper time parametrization (^$t = 0) is equivalent to 
the replacement of the proper time by the physical time t in the kinetic equation 
(8). We extend this requirement to the diffusion (48). In this way we obtain 
the transport equation 

Q+n = --pi'd^n + An,^ + rdf,p^"'wM - o, (49) 

Tc 

where the adjoint of Q^ is in L^{dxdp) . Explicitly, 

-{e + T:^)p-'^{{pwfdP - {wp){wPpP + wPpP)dp + p^vji^wPdp^n. 
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(50) 



If the condition (32) is satisfied then d^^a^'^di, is a non-positive operator as an 
expectation value (33) of the square of a differential operator. It follows (see 
e.g., [39]) that the momentum dependent coefficients a^'' satisfy the positivity 
condition 

af,a^a^"' > 0. (51) 

However, we can prove that the inequality (51) is satisfied under the condition 

e > TTe > (52) 

weaker than (32). We can use the Lorentz invariance of the inequality (51) in 
order to choose a proper Lorentz frame. So, if a^ > we choose a frame such 
that a ~ (oo, 0). In this frame (from cq.(35)) we have 

a^.a^a^'-' = alp-\{e - 7r,)p2 + (e + ^e)(p'w2 - (pw)^), (53) 

where the normalization w'^ = Wq — w^ = 1 has been used. It follows from the 
Schwarz inequality that (53) is non-negative if the condition (52) is satisfied. If 
a^ < then we choose the Lorentz frame such that w = (1, 0). In this frame 

a^.a^.a'"' = (e - n,)p-^(alp^ + pla^ - 2papoao) + 27reP"2(p^a2 - (pa)^). (54) 

(54) is non-negative under the condition (52) again by Schwarz inequality. It 
follows from eqs.(53)-(54) that the conditions (52) are sufficient and necessary 
for a definition of the diffusions (48) and (50) ( the condition (51) is necessary for 
any definition of a diffusion because the diffusion matrix a'^'^ should be equal 
to the correlation function of particle's diffusing paths which is defined by a 
positive probability measure). 

We find the equilibrium solution VIe from the requirement 

{a'^'d^ + rP'"'w^)nE = 0. (55) 

There is a solution (which is the Jiittncr equilibrium [40]), if 

r = /3(e-7r,). (56) 

Then, 

Qe = exp(-/3wp - 7p^), (57) 

where 7 is an arbitrary constant. Eq.(56) can be satisfied for r = if tt^ = e 
which is still in agreement with eq.(52). However, the condition tt^ == e is in 
confiict with the positivity of the measure determining the correlation functions 
of the electromagnetic field ( see eq.(32)). 

If the equilibrium Q,e is known then the form of the friction can be deter- 
mined from the condition of the detailed balance (see an interpretation in terms 
of transition rates in [41] and as Kolmogorov's reversibility condition in [42]) 

f dp^E'^rW = I dpVtE^Wr. (58) 
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This condition is equivalent to the requirement that the generator of the diffu- 
sion is of the form 

g^ = -p^dl + QE'd^a'^'^Ed, (59) 

From eq.(59) it fohows that the friction term of a reversible diffusion must have 
the form 

a^'^(r!^i9^0s)9,. (60) 

Comparing eq.(60) with eq.(48) we can see that the diffusion (48) satisfies the 
detailed balance condition (58)if r is determined by eq.(56). 

6 Discussion and summary 

We have discussed Liouville description of a particle in an electromagnetic field 
coming from a thermal radiation. In such a system there is a distinguished 
reference frame associated with the heat bath. The Poincare invariance of the 
whole system (particle plus the heat bath) can be expressed as a covariance 
with respect to the transformations of the reference frame of an observer. In 
such a description the thermal field becomes a function of the space-time point 
and the reference frame. We have shown that the lowest order term in the ex- 
pansion of the Liouville evolution in powers of the random electromagnetic field 
is described by a diffusion operator which is a generalization of the relativistic 
diffusion of Schay and Dudley. We developed a diffusion theory assuming (af- 
ter Kubo) that the generator of the short time dynamics after multiplication 
by the environmental correlation time determines the diffusion at long time. 
We treat the square of the momentum as a continuous variable. The off shell 
formulation allows to preserve explicit Lorentz covariance. It may have appli- 
cations to quantum field theory (heavy ion collisions) and cosmology (diffusion 
of galactic matter). It seems that the diffusion equation derived from an inter- 
action with a random electromagnetic field could also be inferred on the basis 
of general assumptions of the preservation of the square of the momentum and 
the Poincare invariance of a system described by the phase point {x,p) and a 
fixed vector w which can be interpreted as the velocity of a frame in which 
the particle is studied. Demanding that similarly as in the case of the deter- 
ministic dynamics the probability distribution does not depend on the proper 
time we obtain a diffusive transport equation in the laboratory time. A similar 
transport equation (called a relativistic Ornstein-Uhlenbeck process) has been 
discussed in refs.[31]-[33] without relating it to the physical model of an interac- 
tion with an electromagnetic field. The authors in [31]-[33] investigate diffusion 
of particles in a fluid which supplies a preferred reference frame. However, their 
diffusion does not coincide with the one deflned in this paper. Owing to the 
dependence of a^"^ in eq.(50) on the momenta there is no choice of parameters 
and no choice of the frame which enables to reduce the diffusion (50) to the 
Ornstein-Uhlenbeck diffusion (the Ornstein-Uhlenbeck diffusion appears only in 
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the non-rclativistic limit). It is reassuring that a matheniatical construction of 
the diffusion, discovered first by Schay [17] and Dudley [18] and generalized here 
to realistic systems at finite temperature, is the one which may be realized in 
nature (in particular, in a particle's interaction with a thermal radiation). 
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